##Exerccio 7##

#7.2#
library(MASS)
data(Animals)
anim.m2<-lm(log(brain)~log(body),data=Animals,subset=!(log(Animals$body)>8&log(Animals$brain)<6))
#Ajuste do modelo
anim.m0<-lm(log(brain)~1, data=Animals,subset=!(log(Animals$body)>8&log(Animals$brain)<6))
#Compare dos 2:
anova(anim.m0,anim.m2)

#1. Qual a relao do comando 'anova' acima com:
anova(anim.m2)
# O primeiro faz a ANOVA entre os dois modelos informados, enquanto o segundo compara o modelo informado (anim.2) com o modelo nulo. Ambos geram as mesmas informaes.

#2. Qual a relao entre os valores obtidos por estes comandos:
summary(anim.m0)
mean(log(Animals$brain[!(log(Animals$body)>8&log(Animals$brain)<6)]))
sd(log(Animals$brain[!(log(Animals$body)>8&log(Animals$brain)<6)]))
# o comando "mean" fornece a mdia, mesmo valor da estimativa do desvio padro fornecida pelo comando summary
# o comando "sd" fornece o desvio padro, que  o mesmo valor fornecido pelo resduo do desvio padro fornecido pelo summary.


#7.3#



aves<-read.table("aves_cerrado.csv",header=TRUE,sep=";")
str(aves_cerrado)
aves$carcara[is.na(aves$carcara)]<-0
aves$seriema[is.na(aves$seriema)]<-0
aves$urubu[is.na(aves$urubu)] <- 0
aves

#modelo linear
mod.ce<-lm(seriema~carcara,data=aves,subset=aves$fisionomia=="Ce")
mod.cc<-lm(seriema~carcara,data=aves,subset=aves$fisionomia=="CC")
mod.cl<-lm(seriema~carcara,data=aves,subset=aves$fisionomia=="CL")

#grficos
plot(seriema~carcara,data=aves)
abline(mod.cc, col="blue")
abline(mod.ce,col="red")
abline(mod.cl,col="green")
par(mfrow=c(3,1))
plot(seriema~carcara,data=aves,subset=aves$fisionomia=="Ce")
abline(mod.ce,col="red")
plot(seriema~carcara,data=aves,subset=aves$fisionomia=="CC")
abline(mod.cc, col="blue")
plot(seriema~carcara,data=aves,subset=aves$fisionomia=="CL")
abline(mod.cl,col="green")
par(mfrow=(1,1))

#Coeficientes
coef.ce<-coef(mod.ce) # 7.3921076  -0.1458764 
coef.cc <-coef(mod.cc) # 5.67584746 -0.08262712 
coef.cl<-coef(mod.cl) #7.5833333  -0.4507576 

#anova
anova(mod.ce)
anova(mod.cc)
anova(mod.cl)

p.ce<-0.2401
p.cc<-0.6454
p.cl<-0.1045

# 1.H relao entre o nmero de avistamentos de seriemas e carcars em alguma das fisionomias amostradas?
 No h relao significativa entre o avistamento de carcars e de seriemas em nenhuma das regies do cerrado.

# 2.H diferenas na relao entre o nmero de avistamentos de seriemas e carcars entre as fisionomias do Cerrado?
 H uma relao maior na rea cl do que nas reas cc e ce, pois possui um maior valor de P.